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Abstract 

A finite-dimensional matrix model for the nucleon-nucleon cross section opera- 
tor is used to calculate the dispersive correction to nucleon-nucleus total cross 
sections, and the leading terms in its expansion in the number of inelastic tran- 
sitions in the high-energy limit where the longitudnal momentum transfers can 
be ignored. Results for matrices of different dimensions, but giving the same 
second and third cross section moments, are compared for scattering from ^^^Ph. 
The leading second order terms in the expansion are accurate to better than 
10% , but always larger than the exact result. The 3rd order terms overcorrect, 
but are nearly cancelled by the 4th order terms. The 5th and 6th order terms 
are smaller, and also tend to cancel, but not to the degree of the 3rd and 4th 
order terms. The series converges more rapidly for smaller matrix dimensions, 
and in each order the magnitudes of the terms increase with the dimension of 
the matrix. 



I Introduction 



One consequence of the composite structure of nucleons is a decrease in 
nucleon-nucleus total cross sections due to transitions between different internal 
states of the projectile nuclcon. This decrease can be calulated |^ using an 
operator to represent the nucleon-nucleon total cross section, with the matrix 
elements representing the probability amplitudes for forward scattering transi- 
tions between different states of the nucleon. 

If the longitudnal momentum transfer due to the different masses of these 
states cannot be ignored this dispersive correction can only be calculated by 
expanding in terms of the number of inelastic transitions between states of dif- 
ferent mass. Because of the difficulty of this calculation, and uncertainties in 
the nature of the cross section operator, only the leading term in this expansion 
has been evaluated and then only approximately. (A closely related effect 
in quasi-elastic electron scattering can complicate the analysis of color trans- 
parency effects in these reactions Q.) In this paper we take a first step toward 
testing the convergence of this expansion by using finite matrix models for the 
cross section operator and taking the high energy limit in which the longitudnal 
momentum transfer vanishes. 

Section II reviews the formulas and notation expressing the dispersive correc- 
tions to the hadron-nucleus total crosss sections in terms of an integral over the 
generating function for the cross section probability distribution function. The 
transition expansion is then reviewed in Sect. Ill, where it is shown that each 
term in the expansion can be represented as a sum over products of transition 
amplitudes weighted by a function of the differences among cross sections in the 
different nucleonic states. In Sect. IV the cross section operator is represented 
by a finite dimensional matrix depending on two parameters, one of which fixes 
the second moment of the operator, the other the third moment. With this 
simple form for the cross section operator one can find relatively simple ana- 
lytic expressions for the low order terms in the transition expansion. These are 
evaluated and the results presented in Sect. V. The results are summarized and 
discussed in the concluding Sect. VI. 

II Hadron-Nucleus Cross Sections 

In this paper it is assumed that the energy of the incident hadron is high 
enough that the longitudnal momentum transfers associated with inelastic scat- 
tering can be ignored. Nucleon overlap in the target nucleus is also ignored, 
along with nuclear correlations. For heavy nuclei the total cross section for the 
scattering of an incident hadron from a nucleus of nucleon number A is then 
given by 

C7(A)=C7«(A)-a^(A), (1) 

where 
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aG.-D(A) = 2 J dHG^'^it{A, b)), (2) 
with the Glauber contribution 

= 1 - exp{-t) (3) 
and the "dispersive" correction due to diffractive excitation 

G° =< l\exp{-xt)\l > -exp{-t). (4) 

Here the expectation vahie is in the first, or ground, mass eigenstate of the 
projectile. The dimensionless thickness function t at impact parameter b is 
given by 

t{A,b) = AaT{A,b)/2, (5) 

where a is the projectile-nucleon total cross section and T{A, b) is the usual 
thickness function, i.e. the integral of the nuclear density, normalized to unity, 
along a straight line at constant impact parameter. The dimensionless thickness 
function decreases from a maximum value of approximately {^u / ^-nrQ)A^/^ at 
zero impact parameter to zero for impact parameters well outside the nuclear 
radius R « roA^/^, with ro ~ 1.14/m. Below the values A — 208 and a = 
39.8m6 are used, giving the maximum value of t as about 4.3. The detailed 
shape of t as a function of b depends on the shape of the nuclear density: below 
a Woods-Saxon form will be assumed. 
In Eqn. 4 the operator 

x = a/ <l\a\l> (6) 

is the dimensionless cross section operator. The matrix elements of the cross 
section operator a itself give the cross sections for transitions among the various 
mass eigenstates of the projectile. Thus 

cr =<1|ct|1> (7) 

is the total cross section for the projectile to interact with a single nucleon, 
while , assuming all forward amplitudes are pure imaginary, the total forward 
cross section for diffraction dissociation is 

da/dt\t=o = TT^I < j|<T/47r|l > |2 (8) 

= a^[< l\x^\l > -l]/(167r). (9) 

Since there is diffraction dissociation the operators a and x are clearly not 
diagonal in the space of mass eigenstates. They can, however, be diagonalized 
using their eigenstates \a) , where 
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x\a) — Xa\a). 



(10) 



Then the generating function for the probabihty distribution, needed to evaluate 
Eqn. 4, is 

< l\exp{—xt)\l >= Pq, ea;p(— x^i), (11) 

a 

where 

Pa - I < l|a)P (12) 

is the probabihty of finding the projectile in the scattering eigenstate \a) when 
it is in the mass ground state |1 >. Of course x will have in general a continuous 
as well as a discrete spectrum, so that the sum in Eqn. 11 should be interpreted 
as a sum plus an integral. 

In this paper, however, x will be represented by a finite dimensional matrix 
determined by two parameters, chosen to conform to the following constraints: 

1. The off-diagonal elements < i\x\j > should decrease as \i — j\ increases . 

2. The diagonal elements < j\x\j > should increase, or at least not decrease, 
as j increases. 

Both of these requirements are suggested by non-relativistic wave function mod- 
els: the overlap of the two wavefunctions should decrease with increasing qun- 
tum number separation, and the spatial extent of the wavefunctions should 
increase with increasing quantum numbers. (Related work using finite matrices 
to describe inelastic scattering can be found in |[ 0.) 

Although the longitudnal momentum transfer is ignored in the high energy 
limit used here, it would be reasonably easy to include it along the lines of 
[| if the masses of the excited states were known. This is in fact one of the 
motivations for developing a finite matrix model. 

Ill Transition Expansion 

The operator x can be separated into two components, one diagonal and the 
other off-diagonal (inelastic) in the mass eigenstates: 

x^Xd + xi, (13) 

where 

< i\xd\j >= Sij < i\x\i > . (14) 
Then, just as in the case of time-dependent perturbation theory, 

U{t) = exp{-xt) (15) 



4 



can be expanded in powers of xi : 

oo 

U{t)=e-^^ni + Y,U^/'\t)], (16) 

n=l 

where the nth term is given by the t-ordered integral 

= (-1)" / dtn . . .dtiOitn - tn-i) . . .0{t2 ~ h)xi{tn) . . .xiih), (17) 



with 

xi{t) = e^-'^xie-^"*. (18) 

The diffractive contribution to G{t) can then be expanded in powers of £/, which 
is equivalent to an expansion in the number of transitions between different mass 
eigenstates: 

oo 

G^(t) = ^G^(")(t), (19) 

where 

G^(")(t) = e^* < l|?7j"^(i)|l > . (20) 
Inserting complete sets of mass eigenstates, this can be written as 



= [(_i)"e-Vn!] ^ < l\xi\jn-i > ■ ■ ■ < Ji\xi\l > 
ji,... j„_i 

/(")(y,„_„...,y,J, (21) 



where 

% = i^n - 1)* (22) 
and the functions /("^ are defined as the ordered integrals 



/•I 



/^"••(yn-l, . . . ,?/l) = n\ dui . . .dUnOiUn ~ Un-l) ■ ■ ■0{U2 - Ui) 

Jo 

g-l/n-l(Wn-Mn-l) g-yi(ti2-til) ^23) 

These symmetric functions are normalized to equal one when all arguments 
vanish, and are monotonically decreasing functions of each argument. 

The leading term in this expansion for is well known and often used 
to estimate , especially when the longitudnal monentum transfers are not 
negligible Q: 

G^^'Ht) = (tV2)e-* < ^\^\J >< > f^'^a^n ~ 1)0, (24) 
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where 

f^^\v) = 2[e-y -l + v)/v\ (25) 

One in principle needs to know all matrix elements of x to evaluate G^*-^-* and 
higher order terms, although in practice the sum over j may converge rapidly. 

If all diagonal matrix elements of x are equal then Xd — \^ all y',jS vanish, 
and the above expressions simplify to 

QD(n) ^ (-t)"e~7n!, (26) 

where ii = x — 1 is the shifted operator. 

IV Matrix Model 

The dimensionless cross section operator x determines (t) and its expan- 
sion in inelastic transitions (assuming that the mass eigenstates are known). As 
mentioned above, x may in general have a discrete and/or a continous spectrum. 
(The parameterizations of [0, for example, are based upon a purely continu- 
ous spectrum for x.) Here it is assumed, however, that this operator can be 
approximately represented by a two-parameter matrix of finite dimension N in 
the space of (discrete) mass eigenstates: 

<i\x\j>^S,,ij-l)d + f\^-^\, (27) 

so that the diagonal elements 1,1 + d,l + 2d, . . . increase linearly with i = j 
if d > while the off-diagonal elements decrease geometrically with |i — j| if 
< / < 1. Since 

< lli^ll >=l + .f + .f' + ... + f ('^-i) = (1 - /2^)/(l - f), (28) 

the parameter / is fixed by < l|i;^|l >, which is in turned fixed by forward 
diffaction dissociation according to Eqn. 9. The second parameter d can then 
be determined from 

< > - < l\x^\l >= (2 + d)[f + 2/4 + . . . + (n - l)f'^^-% (29) 

For any N, f and d can therefore be chosen to match the second and third 
moments of x (the only constraint being that d should be non-negative): all 
higher moments and G^{t) are then completely determined. 

This parameterization of x is rather arbitrary and leads to what should 
be considered as a "toy" model. Its main advantages are that it has enough 
flexibility to match any desired < x^ > and < > (with the constraint noted 
above) and that the lower order terms in the transition expansion are relatively 
easy to evaluate. One can think of many other models with these properties: 
to choose among the possibilities one would require a rather complete realistic 
model for the composite structure of the nucleon. 
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In this model the product of matrix elements in each term in the sum in 
Eqn. (21) is just the parameter / raised to some integer power. This multiplies 
a function Z^"-* which is a simple symmetric analytic function of its n — 1 argu- 
ments. Starting with /(^\ given by Eqn. (25), these functions can be calculated 
from recursion relations and power series: 

/^"+'Hyi,2/2,...,yn) = -(n + l)[/("H2/2,2/3,---,2/n)- 

/^"'(2/i,y3,...,2/n)]/(y2-2/i) 
= 1 - (t/1 + 2/2 + . . . + yn)/{n + 2) + . . . . (30) 



V Evaluation for Nucleon-^"*P6 Total Cross Section 

The formulas above can be used to calculate the dispersive corrections to 
any total cross section and their expansions in terms of the number of inelastic 
transitions for cross section matrices of any dimension. Here only the example 
of nucleons on the heavy nucleus will be considered, taking 39.8 mb for 

the nucleon-nucleon total cross section a and using a Woods-Saxon density with 
radius R = 6.75fin and surface thickness parameter = 2.3/m for the 208p^ 
nucleus . (These parameters give an un-corrected Glauber cross section of 3022 
mb.) Here the values < >= 1.25, as in Ref. 1, and < >= 1.9, which 
is close to the smallest value giving non-negative values of the parameter d 
according to Eqn. 29, are used for the 2nd and 3rd moments. (Ref. 1 takes 
< >= 1.75 based upon an approximate analysis of nucleon-deuteron inclusive 
diffractive scattering.) 

The exact integrands AirbG^ in these models, along with the contributions 
from terms of order n = 2, 3, . . . , 6 in inelastic scattering, are shown in Figs. 1 
and 2 for matrix dimensions TV = 3 and N — 6, repectively. The contributions 
from different n's alternate in sign, as required by Eqn. 21, with the positive 
even terms comparable in magnitude to the preceding negative odd terms. This 
can be understood by considering the lowest order terms in powers of the pa- 
rameter which is never much larger than 0.2. A simple analysis shows that 
the leading terms in the G^^^^ and G^^^' are both of order while the leading 
terms in G^^^^ and G^^^^ are both of order Furthermore, the larger the di- 
mension of the matrix the more terms are included in the sum in Eqn. 21, and 
thus the larger the magnitude of G^^"^. The additional terms are in general of 
higher order in however, so this increase converges for large N. 

The results for the cross section corrections after integration over impact 
parameter are shown in Table 1. These results reflect the behavior of the in- 
tegrands: The n = 2 term is accurate to better than 10%, while the 3rd and 
4th order terms range from about 20 to 40% of the exact result, but cancel to 
better than 4%. The 5th and 6th order terms are smaller, and also of opposite 
sign, but do not cancel to the same degree. The sum of the n = 2 through 



7 



n = 6 terms in the expansion gives results which are accurate to better than 1% 
for N=3, but only to about 5% for N=6. (The continuous distribution models 
of Ref. 1 give somewhat higher values, ranging from 199 to 232 mb, for the 
dispersive correction to the total nucleon-^"*^P6 cross sections.) 

The calculations above have been repeated for larger values of < a;^ > up 
to 2.5. The exact valu(^ for the dispersive correction and the magnitudes of 
terms in the expansion all increase with the dimension of the matrix and de- 
crease with increasing < > (presumably because the diagonal elements of 
the matrix become relatively larger as < x'^ > and therefore d increase). The 
near-cancellation of the 3rd and 4th order terms seems to be general, so that 
including these terms never improves the accuracy significantly. 

VI Conclusion 

For the particular matrix models used above, the lowest order term in the 
transition expansion for the diffractive correction to the total cross section gives 
fairly good accuracy, but is always a bit too large. Including higher order terms 
does not improve the accuracy significantly, and in fact adding only the next- 
to-leading 3rd order term always gives poorer accuracy, with a correction less 
than the exact result. 

It would be interesting to repeat these calculations for other matrix models 
to see how general these results are, and to include the longitudnal momentum 
transfers which appear for finite energies. (Preliminary calculations indicate 
that for constant nuclear densities the formulas above can still be used, except 
that the arguments jjj in Eqn. 21 become complex, with imaginary parts de- 
pending on the nucleon's mass spectrum.) Including the longitudnal momenta, 
for example, might change the phases of the terms in the expansion enough to 
interfere with the cancellation of the 3rd and 4th order terms. 
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FIGURES 



FIG. 1 The integrands AnhG^ for the dispersive corrections to the nucleon- 
"^^^Ph total cross section using a 3-dimensiional matrix model as described in the 
text, with < >= 1.25 and < >= 1.9. The solid curve is the exact result, 
the labelled dashed curves the contributions of order 2 thorugh 6 in inelastic 
transitions. 

FIG. 2 The same as FIG. 1, but using a 6-dimensional matrix model. 



TABLE 



TABLE 1 The dispersive corrections to the nucleon-^^^Pfe total cross sec- 
tions for matrix models of different dimensions, including the exact results, the 
contributions from orders 2 through 6 in inelastic transitions, and the sums of 
these five contributions. 



Matrix 
Dimension 


Dispersive Corrections to Total Cross Sections (mb) 


Exact 


Order in Inelastic Transitions 


2 


3 


4 


5 


6 


Sum 


3 


137.64 


143.41 


-30.37 


29.46 


-8.51 


4.58 


138.56 


4 


143.08 


154.82 


-47.13 


47.75 


-20.80 


11.80 


146.44 


5 


145.18 


158.53 


-53.21 


57.02 


-29.60 


18.37 


151.12 


6 


145.82 


159.58 


-54.95 


60.29 


-33.61 


22.13 


153.43 
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